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ABSTRACT 


We evaluate the bit error probability of coded lightwave systems emploving 
direct detection with OOK modulation, and coherent optical systems with FSK 
modulation and noncoherent detection. Block codes, convolutional codes, and con- 
catenated codes are investigated. For direct detection systems, both hard and soft 


decoding are considered. Only hard decoding is employed in coherent optical sys- 


tems. 
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I. INTRODUCTION 


The benefits of forward error-correction coding for radio communications in an 
additive white Gaussian noise (AWGN) channel are well documented [1-3]. With 
the state-of-the-art technology of Application Specific Integrated Circuits (ASIC), 
both block coding and convolutional coding can be implemented at very high bit 
rates. Block coding can be implemented at Gb/s data rates and convolutional coding 
at hundreds of Mb/s data rates. This opens the possibility for employing FEC in 
lightwave systems to compensate for the loss in the !ink budget or to relax the laser 
linewidth requirement in coherent optical systems [4-5]. 

In this study, we explore the benefit of FEC for direct detection lightwave 
systems employing on-off-keying (OOK) modulation, and coherent optical systems 
employing frequency-shift keying (FSK) with envelope detection corrupted by phase 
noise. Our work was derived explicitly from [19]. 

Of practical importance is the use of concatenated codes. Concatenated codes 
allow the use of two relatively short codes with modest error-correcting capability, 
one as the inner code and the other as the outer code to achieve a large coding gain. 
Short codes with modest error-correcting capability simplify decoder design at very 
high data rates and thus reduce the implementation cost. Concatenated codes also 
allow the use of a short constraint length convolutional code as the inner code with 
soft decoding and a short Reed-Solomon (RS) code with modest error-correcting 
capability as the outer code. This scheme can achieve a large coding gain which 
otherwise can only be achieved with a long constraint length convolutional code or 


a long RS code with large error-correcting capability. 











II. BIT ERROR PROBABILITY OF DIRECT 
DETECTION LIGHTWAVE SYSTEMS 


- 


A. OOK MODULATION 

The direct detection OOK receiver model to be analyzed is given in Fig. 1. The 
photodiode detector has a responsivity R = nq/hy (A/W) where 7 is the quantum 
efficiency which is near unity for p ~ 1 —n diodes, ¢ = 1.6 x 107® C is the electron 
charge, h = 6.626 x 10-™ J S is Planck’s constant, and f is the laser frequency. 
The current produced by the photodiode is amplified by a low noise amplifier and 
integrated by the integrator. With soft decoding, the analog samples from the 
integrator are processed by the soft decoder. With hard decoding, the integrator's 
samples are first detected by a threshold detector to provide hard quantized samples 
to the hard decoder. 

In practice, the shot noise and dark currents produced by the photodiode 
are negligible compared to the thermal noise current generated by the low noise 
amplifier. Therefore, they are ignored in our analysis. The complex envelope of the 
OOK lightwave signal incident upon the photodiode detection during a coded bit 
interval T is 


s(t)= VP be) ¢=0,1; O0<t<T (1) 


ae AMPLIFIER ar 
DETECTOR © DECODER be DATA 


Noise 
Figure 1: Coded OOK direct detection receiver structure. 
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where P is peak power, 0) = \/p/(1+ p), & = y1/( + p) and p is the laser extinc- 


tion ratio. The output of the integrator is 
T T 
= (4) |? 9 
i= Rf |si(t)Pae+ f° n(t)dt (2) 


where n(t) is the thermal noise current with spectral density No (A?/Hz). Substi- 


tuting (1) into (2) and denote 
T 
= t)at 3 
N= | n(t)d (3) 
we obtain 
Y¥;=RPTR+N (4) 


Note that N is a Gaussian random variable with zero mean and variance o? = NoT. 
We assume that the soft decoding is done by a maximum likelihood soft decoder 
using the analog (unquantized) samples Y; in (4). For hard decoding, a bit decision 
at the output of the threshold detector employing a threshold setting a is produced 


with channel transition probability as follows: 


are sPr{Yo> a} + 5Pr {vi Za 


2 
1, f(a-—RPTR\ 1, (RPTR-a 
= 30(Srar) +20 (oar) o 
where Q(z) is give by 
Q(z) = Fe [eas (6) 
V2z Je 


The minimum value of p at the optimum threshold setting agp. = RPT(b? + 63)/2 
is 


RPT (b? — 62) 
ee ( ait 


3G (ssNRvr) (7) 
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where 


T, 
SNR = RP(B #8)\/ Xi, (8) 


is the signal-to-noise ratio at the threshold detector input, r is the code rate. and 


T, is the bit interval. 


B. CODING 
Assuming that the decoding is accomplished by a maximum likelihood decoder. 
then the bit error probability of the decoded OOK signals employing block code 


(n,k) is upper bounded by 
PB, < 3 An; P; (9) 


j=l 


where n; is the number of code words of weight j, and P; is the probability that the 
distance between the received sequence and a weight-j code word is less than the 
distance to an all-zero code word 

For a convolutional code with rate k/n, the bit error probability is upper 
bounded by [1-3] 


1 co 
P, < E bs w;P; (10) 
j=dy 


where w; is the total information weight of all code paths of weight 7, d; is the free 
distance of the code, and P; is the probability that the all-zero path is eliminated 
by a path of weight 7 merging with it. 

With hard decoding, the probability P; in (9) and (10) is given by 
8 ( 4 ra ~ py, j odd 


§=(j+1)/2 


P= . : (11) 
Ge pra - py + 2 € )r' —py,  jeven 


s=y/241 








where p is the channel transition probability in (7). Alternatively P, in (9) can be 
calculated by the following upper bound 
P< > det (8) yap (12) 
jai 7 J 
where ¢ is the error correcting capability of the block code. 

For soft decoding, the decoder selects the code word (or code path) at minimum 
Euclidean distance from the received sequence. Thus, P; in (9)-(10) is the probabil- 
ity that the Euclidean distance between the received sequence Y and weight-j code 
word X; is less than the Euclidean distance to the assumed transmitted all-zero 


code word Xo. Thus 
P; = Pr {\|¥ — Xj? < IY - Xoll*} (13) 


where X; is a vector with 7 components, each of value RPTbj, represented bit ones, 
and the remaining components, each of value RPT 83, represented bit zeros; the vec- 
tor Xo represents the assumed transmitted all-zero sequence with each component 
assuming the value RPTb3. The received vector is Y = Xo+N where N is the zero 
mean Gaussian noise vector, each component has variance ao? = NoT. 


We note that by using Y = Xo + N in (12) we have 


P; 


Pr {IN - (X; - Xo)ll? < NIP} 


Pr {\|NII? + IX; — Xoll? — 2Re {(N, Xj — Xo)} < INI}? 


Pr {2Re {(N, X; — Xo)} > [Xj - Xoll”} 


X;—- Xo i ; 
Pr{ Re{ (W ea > slLX, ~ Xo} (14) 


where (N,X; — Xo) is the inner product of the two vectors N’ and X; — Xo. We 


remark that Re {(N,(X; — Xo)/\|X; — Xoll)} is the sum of zero mean Gaussian 
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random variables, each with variance a? = NoT, and (X; — Xo)/||X; — Xol| is a unit 


length vector. Thus, it is a zero mean Gaussian variable with variance 0? = NoT. 
Hence 
X;—- X, 


1 T 
SRP) /i57 | 


Cs 
- o (R280 ere | 
( 


Q 


Q (5SNR jr ) (15) 


where r is the code rate and SNR is given in (8). 

The above results are applied to various codes to study their effectiveness. 
Figure 2 shows the bit error probability P, versus signal-to-noise ratio SNR for 
the (7, 4) Hamming code and the concatenated (7, 4)/(15, 11) Hamming/Reed- 
Solomon codes with the Hamming code as the inner code. With hard decoding the 
Hamming code provides very small coding gain. On the other hand, the coding gain 
at P, = 10-) is about 2.3 dB for the concatenated codes with soft decoding for the 
Hamming code. The coding gain at P, = 107!5 is almost 3 dB as shown in Fig. 3 
for the concatenated (31, 27)/(27, 23) RS/RS with hard decoding where the (27, 
23) is the shortened outer code. 

Figure 4 shows the performance of a convolutional code with rate 3/4 and 
constraint length y = 2 and v = 5. The shorter constraint length »y = 2 is more 
suitable for high data rates. The coding gain at P, = 10~** is about 1.5.dB for v = 2 
and 2.5 dB for vy = 5. Large coding gain can be obtained with the concatenated rate 


3/4 convolutional/(31, 27) RS codes with the convolutional code as the inner code 


eee, ‘ x 
aie MES 


as in Fig. 5. A coding gain of about 2.8 dB and 4 dB can be achieved at P, = 107° 


with soft decoding for the convolutional code with v = 2 and v = 5, respectively. 
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Figure 2: Coded P, versus SNR for direct detection systems employing (7, i 
Hamming code, (15, 11) Reed-Solomon code, and concatenated (7, 4)/(15, 11 
Hamming/ Reed-Solomon codes. 
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Figure 3: Coded PF, versus SNR for direct detection systems employing (31, 
27)/(27, 23) Reed-Solomon/Reed-Solomon codes. 
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Figure 4: Coded P, versus SNR for direct detection systems employing rate 3/4 
convolutional code with constraint lengths v = 2, 5 
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Figure 5: Coded P, versus SNR for direct detection systems employing concate- 
nated rate 3/4 convolutional code/(31, 27) Reed-Solomon codes. 
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III. BIT ERROR PROBABILITY OF COHERENT 
OPTICAL FSK SYSTEMS WITH NONCOHERENT 
DETECTION 


Heterodyne and homodyne (coherent) lightwave systems are largely affected 
by the combined laser phase noise of the source and local laser oscillators. The 
performance deteriorates as the linewidth bit time 67, increases where 6 (Hz) is 
the combined laser linewidth of the source and local lasers and Tj (s) is the bit 
time [6-10]. Coherently detected signals suffer much more from the phase noise 
than noncoherently detected signals. The phase noise which can be modeled as a 
Brownian motion process causes signal fading and broadens the signal spectrum. 
Both effects can be compensated to a large extent by using a wider predetection 
filter to pass more of the signal power and an averaging postdetection filter to 
provide time diversity and thus alleviating the signal fading effect. Such a technique 
was investigated in [11-12] for heterodyne lightwave systems employing ASK and 
wideband FSK with noncoherent detection, and in [13] for homodyne lightwave 
systems. The use of time diversity is equivalent to the use of simple repetition 
codes. Since simple repetition codes can provide a good improvement as reported in 
[11-12], one would naturally expect an even larger improvement when more powerful 
and efficient codes are used. 

In this part we investigate the performance of heterodyne lightwave systems 
using wideband orthogonal FSK with noncoherent detection and error-correction 
codes. We will show that coding can provide a substantial improvement in the face 
of the phase noise and shot noise modeled as a zero-mean Gaussian process with 


power spectral density (PSD) No/2 (W/Hz). Furthermore, coding can relax the BT, 


13 


requirement a great deai thus enabling the use of inexpensive lasers with a much 
larger linewidth than that required by an uncoded system. For total generality we 
consider both time diversity and coding for the system under investigation as shown 
in Fig. 6. We assume that the frequency spacing | fo — f,| is much larger than 1/T" 
where T” is the integration time and is also the sampling time at the summer outputs. 
The time T” is selected such that the code symbol time T = rT, = MT’ where M is 
a positive integer and r is the code rate. Time diversity of the postdetection filter is 
modeled by summing M independent samples of a coded bit after envelope detection. 
The output samples are fed to the decoder which can make a hard decision or soft 
decision on the received codewords. Equivalently, each branch of the demodulator 
in Fig. 6 can be replaced by a bandpass integrator followed by a square-law envelope 
detector [12]. In the following section we derive the channel transition probability. 
We provide a simple derivation as an alternative to the one in [12]. 

For mathematical convenience we adopt the envelope notation of a real signal. 
Thus, the summer outputs of the signal-present upper branch and the signal-absent 


lower branch in Fig. 6 with integration time T’ are 


Yo = aa + jnsoxl” (16) 

k=1 7 
Yi, = Do |rere + irene? (17) 

where Z; is given by : 
Z, = Ap ginas (18) 


2 J(k~-1)T 


and all the in-phase and quadrature noise samples n-i,, Naiz(t = 0,1;k = 1,2,---,M) 
are independent and identically distributed (i.id.) zero-mean Gaussian random 
variables (R.V.) with variance No7'/4. The laser phase noise process 6(t) can be 
characterized as a zero-mean Gaussian process with variance 27 ft [16], where @ is 


the combined linewidth of the source and local laser oscillators. 
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Figure 6: Coded coherent optical FSK receiver structure. 


Since Y; is the sum of squares of 2M i.i.d. Gaussian R.V.s with zero mean and 
variance o? = NoI"/4, Y, is chi-square distributed with the following probability 


density function (pdf) [14] 
fn = aT M-le-nfe’ yy, 20 (19) 


On the other hand, Yo is the sum of squares of M independent Gaussian R.V.s with 
variance a? and mean Re{X;,} and M independent Gaussian R.V.s with the same 


variance o? but with mean J,,{X,}. Here Re{X,} and I,{X,} denote the real and 


15 








imaginary parts of X;, respectively. Thus, Yo is non-central chi-square distributed 


with the following conditional pdf [14] 


l (M~1)/2 /v 
fylyolv) = => (*) en (wote)/207 7 (4) » yr (20) 


207 \v 
where Ij4.1(-) is the (M —1)th-order modified Bessel function of the first kind, and 


v is the vaiue assumed by the random variable V defined as 


M 
V=>-|1%/ (21) 


k=l 


The conditional channel transition probability of error is given by 
p(v) = Pr {¥; 2 Yolv} 
Co foe] 
= [" [ ful fr(olv)dvidyo (22) 


The inner integral can be evaluated as follows: 


00 M-1 k 1 
A fry(yr)dy, = e729? > (3) Ki (23) 


k=0 


Substituting (20) and (23) into (21), we obtain 


e7w M-1 1 00 2 : 
pw) = Sana DB, eu avor de (24) 


where w is the value assumed by the random variable W as defined as 
V 1 
W=-~=-—)(Z,/? (24) 


Using the identity [15] 


roy : : 
[ z/2+n ef 1 (2a,/z )dz = nlavet Ke-n-v=1 LY (-=) (26) 
where 
v ‘ m * ai 
L¥(s) = ZY) ( ee i ml! (27) 





in (24) we obtain 





iis : ie (w/2)™ 
Ae) QM” k=0 2. 3 a (*% jeer m! 
-w/2 M=1 yn Mz os 
= ae Gee | = (28) 
2 n=0 k=n nm 2 
Define 
M 
SNR = A?T,/2No , y= > X;(re/M) 
k=1 
where r is the code rate, c = 2787, and X;, is given by 
ko, a 
X;(re/M) = I vee (29) 
-1 


with w(t) as a zero mean Gaussian process of variance t; we can express the condi- 


tional channel transition probability in (28) as follows: 


e~(SNR)ry/2M M-1 [((SNR)ry/2M]" M-1 k+M-1 1 
The channel transition probability is given by 
p= EXply)} = | rly) fr(u)dy (31) 


where fy(y) is the probability density function of y and can be calculated by the 
method given in (12]. Figure 7 shows fy(y) for M = 2,c = 1, and c = 4 for the 
uncoded case. 

The transition probability p in (31) is applied to (9)-(12) to evaluate the 
effectiveness of various codes for coherent optical systems. Figure 8 shows the bit 
error probability FP, versus the signal-to-noise ratio SNR for M = 1 and-c = 1 (small 
phase noise) with hard decoding (15, 11) and (31, 27) RS codes and (31, 27)/(27, 23) 
concatenated RS/RS codes. It is seen that coding provides an excellent performance 


as compared to an uncoded system, especially for the concatenated code. When 
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M = 2, that is, the bandwidth of the integrators is twice the coded rate. more 
signal power passes through the demodulator and the performance of an uncoded 
system as shown in Fig. 9 improves remarkably over an uncoded system for M = 1 
in Fig. 8. In this situation the improvement in perfor::.ance when RS codes are 
employed is less than that in Fig. 8 but the improvement is still considerable. 
especially for the concatenated (31, 27)/(27, 23) RS/RS codes. When the phase 
noise is large as shown in Fig. 10 for c = 4, the use of both RS codes and larger 
integrator bandwidth M = 2 achieves considerable improvement for the (15, 11) and 
(31, 27) RS codes and the concatenated (31, 27)/(27, 23) RS/RS codes. Figures 11 
and 12 show the performance of rate 3/4 convolutional codes with hard decoding for 
M = 2, andc=1 and c = 4, respectively. It is observed that when the constraint 
length is small, v = 2, the improvement is small. The performance of rate 3/4. 


v = 5 convolutional codes is slightly less than that of RS codes (15, 11) or (31, 27). 
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Figure 7: Probability density function of y for M = 2,c=1, andc=4. 
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Figure 8: Coded P, versus SNR for coherent optical systems employing the (15, 
11) and (31, 27) Reed-Solomon codes and the concatenated (31, 27)/(27, 23) Reed- 
Solomon /Reed-Solomon codes with M = 1 and c=1. 
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Figure 9: Coded P, versus SNR for coherent optical systems employing the (15, 
oA and (31, 27) Reed-Solomon codes and the concatenated (31, 27)/(27, 23) Reed- 
Solomon /Reed-Solomon codes with M = 2 and c= 1. 
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Figure 10: Coded P, versus SNR for coherent optical systems employing the (15, 
A and (31, 27) Reed-Solomon codes and the concatenated (31, 27)/(27, 23) Reed- 
Solomon/Reed-Solomon codes with M = 2 and c = 4. 


22 





Probability of Bit Error 


10 3611 12 3 #14 #2415 6% #17 ~«2«18~«49 
SNR (dB) 


Figure 11: Coded P, versus SNR for coherent optical systems employing hard 
er pay convolutional code with rate 3/4 and constraint length 2, M = 2,c = 1 
and c= 4. 
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Figure 12: Code P, versus SNR for coherent optical systems employing hard 
decoding convolutional codes with rate 3/4 and constraint 5, M = 2, c = 1 and 
c= 4, 





IV. CONCLUSIONS 


From our investigaion we conclude that error correction coding can improve 
the performance of both direct detection and coherent optical systems. The im- 
provement is considerable when concatenated codes are used. From the practical 
point of view, short RS codes with modest error correcting capability can be con- 
catenated to provide large coding gain and at the same time can be implemented at 
high data rates. Convolutional codes with soft decision also provides good coding 
gain at moderate data rates for direct detection systems. It remains to be seen 
how much coding gain can be obtained with soft decision convolutional codes for 
coherent optical systems. 

Our contributions in this thesis are the numerical results obtained from the 


analytical work in [19] as shown in equations (1-31). 
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APPENDIX 


The following MATHEMATICA (Stephen Wolfram, 1993) and MATLAB (The Math 


Works, 1993) codes offer the formulations for: 
The convol: .ion method for the determination of the probability density func- 
tion 


The computation of bit error probability for direct detection and coherent 


detection 


1. Obtaining of pdf of y by M-fold convolution 


£18 [ies cores Td) 0 cteeea eek eee een aS Convolution method for finding pdf of y 
eek e oe ee eee Input of pdf of x 

Z2Q™X1i0 | AA Sash sd Mog eerste d elon t eo ta wer xi,x2 are i.i.d 

x320:0.02:2; 

yi=0.02*conv(x1,x2); we eee eee eee 0.02 is integer as interval 


EE er eee 2nd convolution 


2. Bit error probability for direct detection 


2.1 Uncoding case 


function y=q(x) ke eee eee Definition of q(x) in MATLAB 
y20.S*erfc(x/sqrt(2)) «wwe eee eee Use error function defaulted in 
MATLAB 


snrdb#6.5:0.02:12.1; 
snr=10.*(snrdb*0.1); 
pusq(snr*0.5); he ee ee eee Uncode case 


2.2 Concatenated convolutional hard and soft/RS coding 


Mi=input (’ input inner convolution code order as n,k,v,df,wd,wdt+1...’); 
eee ee ree Need to input : 
two matrix,one for RS code(outter), 
one for convolution (inner) 
M2=input(’input outter RS code as n,k’); 
n1=M1(1); ki#M1(2); n2=M2(1); k2=M2(2); 
rizki/ni; r2#k2/n2; 
m2=log2(n2+1); t2=(n2-k2)*0.5; 
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rerier2; 
p=q(snr*r“0.5*0.5); Equation (7) 
; Solve five pj values 
Hard case 
Soft case 
for d=0:4 
pj=0; 
j=Mi (4) +d; 
if rem(j,2)==1 Check odd or even case 
for i=(j+1)/2:j 
pipjbino().,i)9p- "2. *Cd-p) -"0973)3 
en 


else 
for i=j/2+1:j 
pj=pj+bino(j,i)*p.*i.*(1-p).*(j-i); 
end 


pj=pj+0.5*bino(j ,0.5*j)*p.°(0.5*j).*(1-p) .~(0.5*#j); 
d 


en 
pji=(pji;pjl; Hard case’s pj 
ppjjzq( (r*j)*0.5*snr#0.5); Soft case’s pj 
pj2=(pj2;ppjjJ; 
end 
pb_conv=0; Solve for pb for hard case 
ppbb_conv=0; For soft case 
for h#1:5 5 is from wd to wdt4 
pb_conv=pb_conv+M1i(h+4)*pji(b,:); 
ppbb_conv=ppbb_conv+M1 (h+4) *pj2(h,:); 
end 
pb_conv=pb_conv/k1; 
ppbb_conv=ppbb_conv/k1; 
ps*1-(1-pb_conv) . “m2; Solve for ps 
ppss#i-(1i-ppbb_conv) . “m2; 
pb=0; ppbb=0; Solve for pb 
for jxt2+i:n2 
pb=pb+(j+t2)/n2*bino(n2,j)*ps.*j.*(1-ps) .~(n2~j); 
ppbb=ppbb+(j+t2) /n2*bino(n2,j)*ppss.~j.*(1-ppss) .*(n2-j); 
end 


pb=pb*(n2+1) /(2*n2) ; 
ppbb=ppbb* (n2+1) /(2#n2) ; 
ss*linspace(1,len,51); Plot statement 
semilogy (snrdb(ss), pu(ss),’*’,snrdb,pb,’--’ ,snrdb, ppbb) 
axis(({snrdb(1) snrdb(len) 10°(-15) 1]) 
3. Bit error probability for coherent detection 


3.1 Uncoding case 


Input pdf value of y 
Wanted SNR 








yi=Table[y,{y,0.01761,1.9,0.02}]; 

pi=Table(Sum[(((S#y)/(2#M))~ n)/(n!) 

*Sum [Binomial [M+k-1,k-n] 

#(27(-k)) {k,n ,M-1}] ,{n,0,M-1}] 

*Exp[-(S*y/(2*M))]/(2°M), 

{y,0.01761,1.9,0.02}]; 

p2=Table[{y1[Ci]] ,pillil]*y11[[i]]}, 

{i,1,Length[y1]}]; 

p3=Interpolation[p2, InterpolationOrder -> 1] 
ere eT ee ee Connect each point 

p4=NIntegrate[p3[y] , {y,0.01761,1.89761}] 


3.2 concatenated RS/RS coding 


yliet colcacose. PS. ee ares hee For coding case’s pdf of y value 
For(M#1; g21; ni#31; n2*27; ki=27; k2=23; rci=k1/n1; 

re2=k2/n2; SI*10§ 2 | dh Fue eHs Meek ts Inputs of SNR and code rate 
S#Si*rci*re2, S +35, 

yi=Table[y, {y,0. 04283 ,0.97772,0.02}) ; 

pi=Table[Sum[(((S#y)/(2#M))~ n)/(n!) 

*Sum (Binomial [M+k-1,k-n] 

#(27(-k)) ,{k,n,M-1}] ,{n,0,M-1}] 

*Exp[-(S#y/(2#M))1/(2°M), 

{y,0.04283,0.97772,0.02}] ; 

pa-Table({yiC(il] piiLil]*y11([i]]}, 

{i,1,Length[y1}}]; 

p3=Interpolation[p2, InterpolationOrder -> 1]; 

p4*NIntegrate([p3[y] , fy ,0.04283,0.96283}] ; 

Print [p4J] 


3.3 Convolutiona) coding 


p=( 

len=length(p) ; 

M=#(4,3,2,3,15,104,540,2520,11048);  .......... Matrix order as 
n,k,v,df,wd,wdt1,wd+2,... 


TOMCQ) (MCLs. cick ee ae eee ee ee es Solve for convolution code in 
matrix 

DIS oy cee hunts tos Te ipivenetatews Solve five Pj values 

For G80: occ cece cere reer e cn eceees d i.e., df which value depend on 

weight wd+4 

pj=0; 
pj2M(4)+d; 

Sf vemQ,.2)881 i bee discs cnwdincaas Check odd or even case 

for i=(j+1)/2:j 

pj*pj+bino(j,i)*p.“i.*(i-p) .*(j-i); 

end 
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else 
for i=j/2+1:j 
pj=pj+bino(j,i)*p.~i.*(1-p) .“(j-i); 
end 
BI Pato pep ane ts s0.88) ep k02b29)-* Cop) 2 “(0-52))) 
en 
pji=(pji;pjl; 
en 


pb=pb+M(h+4) *pj1(h, :); 
end 


pb=pb/M(2) ; 
PLOCKE;BD) ° 2 ls SU See ehiceee sees Plot statement 
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